Virialization-induced curvature as a physical explanation for dark energy 
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The geometry of the dark energy and cold dark matter dominated cosmological model (ACDM) is 
commonly assumed to be given by a Friedmann-Lemaitre-Robertson-Walker (FLRW) metric, i.e. 
it assumes homogeneity in the comoving spatial section. The FLRW approximation is expected to 
fail at (i) small distance scales and (ii) recent epochs. We use the virialization fraction to quantify 
(i) and (ii), which approximately coincide with each other on the observational past light cone. For 
increasing time, the virialization fraction increases above 10% at about the same redshift {r^ 1-3) 
at which 57a grows above 10% (« 1.8). Thus, instead of non-zero Q.a, we propose an approximate, 
general-relativistic correction to the matter-dominated (fim = 1, Ha = 0), homogeneous metric 
(Einstein-de Sitter, EdS). A low-redshift effective matter-density parameter of fij^f (0) — 0.26 ±0.05 
is inferred. Over redshifts < z < 3, the distance modulus of the virialization-corrected EdS 
model approximately matches the ACDM distance modulus. This rough approximation assumes 
"old physics" (general relativity), not "new physics". Thus, pending more detailed calculations, we 
strengthen the claim that "dark energy" should be considered as an artefact of emerging average 
curvature in the void-dominated Universe, via a novel approach that quantifies the relation between 
virialization and average curvature evolution. 

PACS numbers: 98.80.Jk, 04.20.Gz, 02.40.-k 



I. INTRODUCTION 

It is widely believed that a valid general-relativistic in- 
terpretation of recent extragalactic observations, in par- 
ticular faint galaxy number counts [e.g. [H-Qj gravita- 
tional lensing [e.g. supernovae type la magnitude- 
redshift relations [e.g. @), and the Wilkinson Mi- 
crovifave Anisotropy Probe (WMAP) observations of the 
cosmic microwave background 7], is that the present- 
day Universe is dominated by a non-zero dark energy 
term, modeled in the simplest case by a cosmological 
constant with today's value r^Ao ~ 0.68 [S;]- This in- 
terpretation is a consequence of forcing the Friedmann- 
Lemaitre-Robertson- Walker (FLRW) model [ol-ITsj. an 
exact, homogeneous solution of Einstein's equations, onto 
the observational data. What are the key assumptions in 
applying the FLRW model to the data? 

One key assumption is an applied mathematics hy- 
pothesis, called the "Cosmological Principle" that we 
rephrase here in a weaker form than usual: a synchronous 
space-time foliation of the Universe should exist accord- 
ing to which each spatial section can be approximated 
by a constant-curvature metric on some assumed large 
scale of homogeneity, and the evolution of this metric 
is given by the homogeneous-isotropic FLRW solution. 
Since the real Universe is (obviously) inhomogeneous, a 
second assumption is required: that the cosmic web of 
filaments, clusters of galaxies and galaxies themselves in- 
duce only small perturbations of the perfectly homoge- 
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neous geometry of the FLRW background. We may think 
of the background as a template space-time, the validity 
of which should be questioned. (For the question of which 
background to use and the construction of template met- 
rics employed for the interpretation of observations, see 

iMiieiJ-) 

At what times and length scales are these two as- 
sumptions most questionable? They are least accurate 
at recent times (redshifts z ^ 3) — since galaxies and 
large-scale structure have mostly formed recently — and 
at small (<§; c/TJq = 3 Gpc) length scales — where 
structures have had sufficient time to become non-linear 
and to be bound as virialized objects. The assumption 
that these structures imply small perturbations of the 
geometry of the background may be valid for metric per- 
turbations, but it is violated for the derivatives of the 
metric, in particular for the intrinsic curvature [17i] . 

On the observational past light cone, the recent epoch 
and small length scale regimes, at which we expect that 
the FLRW solution may fail, approximately coincide. 

In the FLRW model, Vl\ varies with time. At what 
epochs is J^a significantly non-zero, assuming that the 
present-day value is JIao — 0.68? Dark energy is signifi- 
cantly non-zero at recent epochs (z ^ 3). 

Hence, the FLRW model should reasonably be ex- 
pected to be a bad approximation on the same scales 
at which J7a inferred from applying the FLRW model 
is significantly non-zero. Since general relativity is well- 
established empirically on many different scales, the con- 
servative scientific approach is to assume that non-zero 
dark energy is an artefact of trying to apply the FLRW 
model in domains where it is expected to fail, unless or 
until proven otherwise. 
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FIG. 1. (colour online) Redshift evolution of Q,a [Eq. H}] 
in the ACDM model compared to that of the fraction of viri- 
alized non-relativistic (non-baryonic plus baryonic) matter, 
/vir(z), in EdS 256^-particle Virgo Consortium [l l, [l9 [ A^- 
body simulations of box sizes as labelled (see Sect. Ill For 
comparison with Qa, 2/vir is also shown for the 240 Mpc 
simulation. 

Can the relation between the scales of the expected 
failure of the FLRW model and the significantly non-zero 
values of FLRW-inferred JIa be quantified? An obvious 
statistic to quantify the inhomogeneity of the Universe is 
the fraction of non-relativistic matter (baryonic and non- 
baryonic matter together) that is contained in virialized 
(gravitationally bound) systems at a given epoch, fvh-iz). 
In a fiat FLRW model, ignoring radiation density. 

Figure [T] shovifs that on a linear redshift scale, the 
epoch of virialization in an Einstein-de Sitter (EdS, 
f^mo = Ijf^Ao ~ 0) model with a cold dark mat- 
ter (CDM) initial power spectrum (EdS-CDM) coin- 
cides closely with the epoch at which JIa becomes non- 
negligible in the Concordance Model ([20|, ACDM), with 
the details depending on the limitations of the TV-body 
simulations (box size, mass resolution, group-finding al- 
gorithms, etc.). The EdS-CDM virialization fraction and 
the FLRW-inferred dark energy f^A evolve similarly, to 
much better than an order of magnitude. Thus, the 
FLRW-inferred J^a dominates the Universe in a way that 
approximately imitates the degree to which the FLRW 
geometry should fail, shown in Fig. [1] for an EdS-CDM 
model. 

How can the virialization fraction and the FLRW- 
inferred dark energy parameter be more precisely re- 
lated? Since virialized matter occupies about 1/100-th 
to 1/200-th of the volume that it would occupy if dis- 
tributed according to the FLRW background density (e.g. 



App. A in pi!]), the volume- aver age of the 3-spatial cur- 
vature should become dominated by the curvature corre- 
sponding to the volume occupied by the remaining non- 
virialized matter, generating, on average, negative spatial 
curvature (if the initial spatial curvature is zero) . 

In an EdS model the virialization overdensity is, ac- 
cording to a top-hat estimate using the scalar virial the- 
orem for isolated systems, 

5^;, = IStt^ « 178 , (2) 

and for low density, zero dark energy models it is bounded 
below by Stt^ [21[. Since the curvature is effectively neg- 
ative, the density of the non-virialized matter inferred 
from assuming zero curvature is an overestimate. Since 
matter mostly fiows out of voids, the rate of expansion 
(Hubble parameter) in the majority of the volume, i.e. 
in non-virialized low-density regions, is higher than the 
background model average expansion rate. Again, the 
matter density parameter, inferred from assuming a ho- 
mogeneous Hubble parameter, is overestimated in com- 
parison to calculating it in proportion to the critical den- 
sity of the non-virialized region itself. Thus, the degree 
of negative curvature is underestimated both due to cur- 
vature itself and due to deviations in the expansion rate, 
which can be condensed into a single A'-matter parameter 
(see below). 

The scalar averaging approach to relativistic, inho- 
mogeneous models of the Universe [e.g. [13, [22| - [25| . 
focuses on the kinematical backreaction, ilg(z) = 
— Qv/idH^), and the average curvature parameter 
ri^(z) = —{TZ)d/{6H.^), where Hj) is the volume- 
averaged Hubble rate within a spatial compact domain 
v. The former is an averaged expression of extrinsic cur- 
vature invariants, which can be interpreted kinematically 
and depend (in general) on the variance in the expansion 
rate and the averaged rates of shear and vorticity, while 
the latter represents the average scalar curvature (3-Ricci 
curvature) over the domain V. The overall backreaction 
effect can be condensed into an A'-matter parameter, de- 
fined i^'^{z) :— riQ{z) + ilJj^{z) and obeying the Hamilton 
constraint on the averaged density parameters 

nZ{z) + n^iz) + n^{z) = i, (3) 

where fl^ := A/{3H^) is the volume-averaged dark en- 
ergy parameter. 

Earlier calculations, starting with a homogeneous EdS 
background model, suggest that the peculiar curvature 
parameter, i.e. the deviation of the total averaged 3-Ricci 
curvature from a constant-curvature model, (z) := 
-Wv/6H^, with Wv := - 6fc/a^ has a much 

stronger (a factor of about 5) effect than the kinemat- 
ical backreaction r2g(z) (cf. Figs 3, 4 in [22|). Thus, the 
use of virialization to calculate more realistic estimates 
of the average scalar curvature and its evolution should 
provide an approximate, general-relativistic space-time 
model that is more accurate than the FLRW model. This 
effective-metric model is not expected to satisfy the Ein- 
stein equation on any given (recent) time slice [16j . Thus, 
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the effective metric presented below ((38)) is unlikely to 
be consistent with the Lemaitre-Tolman-Bondi (LTB) 
model p6l - l28j . which has recently been parametrized 
against observations [1^ [13] ■ 

The approximation method proposed here, for an EdS 
background model, is described in Sect. HIl and the results 
are presented in Sect. IIIII An OCTAVE script for making 
the calculations and plots is provided in the source pack- 
age for the preprint version of this paper [sij. Conclusions 
are given in Sect. [V] 



II. THE VIRIALIZATION APPROXIMATION 
A. A'^-body simulations 

"Dark" (i.e. baryonic plus non-baryonic) matter halo 
merger history trees have been calculated since 1992 from 
iV-body simulations [s^lssj and by semi- analytical meth- 



ods [3J|. Estimates of the virialization fraction are im- 
plied by these calculations, e.g. Table 6 in [35^. Here, 
the A^-body simulations used are Virg o Consortium 256^- 
particle (3-torus) simulations [ll H^, Hlj for an EdS- 
CDM model with h = 0.5, and normalization in the 
mean density fluctuation as ~ 0.51 at Sh"^ Mpc, where 
the Hubble constant, i.e. the zero redshift Hubble pa- 
rameter, is written Hq = lOO/i km/s/Mpc. One simu- 
lation has comoving box size 239.5 Mpc and mass 
per dark matter particle (implicitly, baryonic and non- 
baryonic together) of to = 2.27 x 10^^h~^MQ, the other 
84.55 h-^ Mpc and m = 1.0 x IO^^H-'^Mq. The length 
scale of the latter simulation is small, so the simulation 
is only used for comparison in Fig. [T] A friends-of-friend 
group finder fof-1.i[3^ at 0.2 times the mean interpar- 
ticle separation, for a minimum of 8 particles per group, 
i.e. 1.8 X 10^^h~^MQ, is used independently at each out- 
put redshift z to detect virialized objects at that red- 
shift, giving the number of virialized particles TVvir and 
the complement, i.e. the number of non-virialized par- 
ticles N — TVvir, where N — 256'^ is the total number of 
particles. A different group finder, such as a spherical 
overdensity (SO) group finder, wou ld give somewhat dif- 
ferent results to those found here '37H39l|. but since we 
are interested in the total virialized mass, the dilemma 
of whether to define a slightly overlapping pair of haloes 
as a single halo (FOF) or a pair of dynamically distinct 
haloes (SO) is insignificant for the present work. We de- 
fine the virialization fraction 



iVvir(g) 

N 



(4) 



We spline interpolate between simulation output times. 



B. Effective cosmological parameters 

The effective, i.e. volume-averaged, matter density pa- 
rameter combines the matter density parameter in the 



non-virialized region with that in the (at late times) 
much tinier virialized region. If we think of the particles 
in their original comoving positions in the background 
model, then the volumes occupied by the two regions are 
approximately in the ratio (1 — /vir) : /vir- However, 
taking into account the actual situation of a virialized 
region, we notice that the volume occupied by the viri- 
alized matter has shrunk by a factor of about Svir ©j 
so in the homogeneous model (with no local nor comov- 
ing global curvature changes), the volume ratio increases 
to (1 - /vir/i5vir) : /vir/i^vir- FoUowiug, e.g., [H, E^] , let 
us label the non-virialized and virialized regions £ (for 
"empty" ) and A4 (for "massive" ) , respectively, and their 
disjoint union V := A4 U £ . We now establish a general 
formalism in Sect. IIIB H derive formulae for the dark- 
energy-free, stable clustering case in Sect. Ill 1^21 and de- 
scribe how to use these in Sect. IIIB 31 



1. General formulae 

Writing | for the volume of a spatial region J^, let us 
define Xm ■= \M \ / \T>\. This corresponds to the virial- 
ization volume fraction, with A;vi:=/vir/<5vir- As in (25) 
of [23j . we then have the linear combinations 

Hv = XmHm + {1- Xm)H£ (5) 
= Am + (1 - Aa^) (6) 
{n)^ = Xm {n)^ + (1 - A^,) {U), . (7) 

As shown in [23l |. the the volume-averaged Hamiltonian 
constraint, 
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A 



(8) 



(9) 



leads to the kinematical backreaction on P, 

Qv = XmQm + {'^- A>i) Qe 

+6Xm (l - Xm) (Hm - Hsf . 

As in j40| , we may now define cosmological parameters 
on V using Similarly to the FLRW model, these are 
derived by dividing ([8|) by 3H^. For a generic spatial 
domain J-', which may be any of I?, A4, or £, we define 
the cosmological parameters 



8ttG I , 



3H. 



"a 



A 



ml 



ml 



(10) 



The choice of dividing by Hi independently of the choice 
of T is motivated by the stable clustering hypothesis f4lj , 
according to which the averaged expansion rate in M. 
(virialized regions) is zero, i.e. Hm ~ 0. 

With these definitions, the Hamiltonian constraint (|8]) 
on T is 



Hi 
Hi 



(11) 
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Thus, the dimensionless parameters fl defined this way 
sum to unity on the combined domain I?, but on M or 
£ are only constrained to sum to a non-negative value, 
which is determined by comparing the region's squared 
expansion rate to that of the combined domain 2?. 



2. Dark-energy-free, stable clustering case 

Let us assume no dark energy, i.e. A = 0, and the 
stable clustering hypothesis for the M regions, Hj^ ~ 0. 
Thus (ISl) becomes 



of equations is 



- 



(1 — Xm) f^m + ^M^m 
Xjli^X^ + (1 — Xm) 



1 



1- A 



M 



(16) 

(17) 

(18) 
(19) 

(20) 



where (PO)) is redundant, since it is implied by (ITC)) -([n 
Thus, there are four independent algebraic relations with 
six unknowns, given the virialization volume fraction Xm ■ 



H-D — (1 — Xm) Hg 



(12) 



Using (112]), dH), 0, dS]), and (Uni) we find that the average 
characteristics are related: 



{1-Xm)^L 
(1- 



XM^m ; 



= ^M^Q + (1 ^ ^m)^q — - 



Xm) ^ 



£ . 



Xm 



(13) 



The Hamiltonian constraint in the form (|lip gives three 
independent equilibria. 
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(I-Xm)' 



(14) 



Together, (|T3|) and p4| consist of, at a given redshift z, 
six algebraic equations with nine unknowns, where the 
virialization volume fraction Xm i^) is estimated by cal- 
culating /vii (z) from a background model ^ and by us- 
ing the appropriate analytical value of 5vii for that model. 
This set of equations should be used along with the void- 
dominated expansion rate , since H-p is used in defin- 
ing these Q parameters. These formulae are valid at each 
redshift. 

To retain generality while simplifying the algebra, we 
introduce the sums 



(15) 



so that we can consider X matter to represent the devi- 
ation of the averaged parameters from FLRW behavior. 
Thus, six variables are reduced to three, while equations 
([7]) and ^ become a single equation. The resulting set 



3. Background model 



The equations in Sect. IIIB II could, in principle, be 
solved in a background-free way. Here, we use a back- 
ground FLRW model to approximate the void region 
matter density parameter 



(1 - /vir) {py- 



= (1 - /vir) 



(21) 

where (p}^^,f2jjf, and H^^ are the FLRW mean den- 
sity, density parameter, and expansion rate, respectively. 
Here, we have assumed that since /vir < 1 and <5vir ~ 178 
[EdS case, ([2])], we have Xm <C 1, so that the non- 
virialized matter occupies approximately the full volume. 
The virialized matter density parameter can also be es- 
timated from the background model. 



M 



8ttG , , 



SttG 


<5vir (P) 






<5vir (y 





bg 



(22) 



Thus, substituting ([21]) and (|22j) into (HH) gives 

/ijbg 



f)f (z) 



•D 



1- 1^(1 -/vir) 



i/bg 



(23) 



where the 2?-averaged parameter has been relabeled as 
an effective parameter. For percent-level accuracy, the 
first line in ([25)1 is necessary, though not sufficient. 

When assuming a background EdS model, the effec- 
tive expansion rate implied by the stable clustering hy- 
pothesis, i.e. using (jl2p . combines the background, time- 
varying expansion rate H^^{z) with the peculiar expan- 
sion rate expressed as a velocity difference normalized by 
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a comoving spatial separation, i7p°™(z), 

H''''{z):=Hv (24) 
= (1-Aa,)[F'^s(z) + 7J™-(z) a-i] , 

where the a^^ factor converts the peculiar expansion rate 
from comoving to physical length units. At redshifts 
prior to the virialization epoch, i.e. for ^ 3, we must 
have 

i7°^(z) « H^^iz). (25) 
This is given by the homogeneous Friedmann equation. 



AO 



(26) 



for background model zero-redshift parameters includ- 
ing (in general) the Hubble constant (Hq^), and the dark 
energy (ri^o)' curvature (f^^g), and matter density (^^J^o) 
parameters. In a comoving-rigid void model (i.e. a stan- 
dard FLRW model), these background parameters cor- 
respond to zero-redshift parameters. In the virializa- 
tion approximation, these parameters only exist as hy- 
pothetical projections from high redshift to an idealized 
low-redshift state; they are not physically realized. In 
the present work, we only consider an EdS background 
model, so (1^51) simplifies to 



(27) 



In order to match low redshift estimates, the background 
model Hubble constant H^^ should be set so that the 
virialization approximation zero-redshift value 



/vir(O) 



rrcom 
pec 



(0) 



(28) 



[cf. ((M)) and (l?7l) ] matches recent low redshift estimates 
of i?o (Sect.lHSl). 

Estimating Hp°™{z) from background iV-body simula- 
tions would be difficult with standard FLRW simula- 
tions, since the average velocity difference is calculated 
along a straight closed curve 7 across the simulation box, 
i.e. it is zero by construction. Thus, here, we estimate 
ffp°™(0) from observations (Sect. UrCj) . Determining the 
redshift dependence Hp°^{z) directly from observations 
would be more difficult and model-dependent. However, 
we have two known constraints: ([25)1 and ([^5)) must be 
satisfied at the high and low redshift limits, respectively. 
At the high redshift limit, a smooth transition at high 
redshifts seems physically reasonable. Moreover, virial- 
ization cannot occur without reducing the matter density 
in the voids. As the voids become more and more empty, 
the imbalance in gravitational potential between the cen- 
tre and edge of a void becomes stronger and stronger. 
Thus, it is physically likely that i/p°™(z) increases mono- 
tonically as fvir{z) increases (i.e. as z decreases). Here, 



we propose a functional form for Hp°™{z) proportional 
to the virialization fraction, i.e. 



rrcom/ \ rrcom 

pec \ ) pec 



(0) 



fvirjz) 
/vir(O)' 



(29) 



By definition, this formula satisfies the high and low 
redshift limits, has a smooth high-redshift transition, 
and has Hp°™(z) increasing monotonically as fvir{z) in- 
creases. Checking the detailed shape of this function 
should be considered as both an observational and a the- 
oretical test of the hypothesis that the EdS model is the 
correct background model. Neither task is trivial, how- 
ever, since both need to be done in the framework of 
the scalar averaging approach (e.g., using the relativis- 
tic Zel'dovich approximation [2J, |25|), rather than in the 
perturbed FLRW framework. 

Since we assume zero dark energy in our background 
model and we approximate il^{z) as small, p^ . ([3]), and 
(|15|) give the effective sign-reversed curvature parameter 



n^^{z) = i-n'i^{z) 



(30) 



where I?-averaged parameters are again relabeled as ef- 
fective parameters. The effective comoving curvature ra- 
dius at a given epoch can now be written 



1 



aH^^{z) 



(31) 



where for simplicity, we assume that rj^f > Vz. 

Using (P7)) . we can rewrite the first line of ([25)1 for the 
EdS case as 



17f(z) 



/v 



-(l-/vh-) 



(32) 



i.e., an effective matter density parameter that is clearly 
less than unity during the virialization epoch. 



C. Observational assumptions 

Two low-redshift limit observational estimates are 
needed in order to evaluate (|24p . We adopt 



i/'=«(0) = 74.0 ± 1.6 km/s/Mpc 



(33) 



using the recent low-redshift Hubble constant estimates 
of Ho = 73.8 ± 2.4 km/s/Mpc ^42] and Hq = 74.3 ± 
2.1 km/s/Mpc HI]. 

To estimate the present value of the peculiar expansion 
rate i?p°™(0), we divide the typical infall velocity Winfaii 
around rich clusters of galaxies by the typical void radius 
-Dvoid/2, where both are for low redshift data. Out to 
a distance of 130h^^ Mpc and at 10 or more degrees 
above the Galactic Plane, the median size of eight voids 
found in the IRAS Point Source Catalog redshift survey 
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(IRAS/PSCz) is i?void/2 = 28.3 /i"! Mpc [i^. In Data 
Release 5 of the Sloan Digital Sky Survey (SDSS), a much 
larger sample of 222 voids found in the redshift range 
0.04 < z < 0.16 in a solid angle of about 2300 deg^ has 
a mean comoving void radius of 



Z^void/S « 25 ± 2 h-^ Mpc 



(34) 



fp5|: standard error by inspection of Fig. 4). 

Infall velocities are not normally derived from observa- 
tions with the aim of estimating Hp°™. The Cluster and 
Infall Region Nearby Survey (CAIRNS) study of nine 
low-redshift rich clusters gives a caustic outline for infall 
velocities in front and behind the clusters of about 



Winfaii ~ 2cr„, 



(35) 



where ay is the line-of-sight velocity dispersion of clus- 
ter galaxies (Fig. 4, [i^, within one sky-plane virial ra- 
dius of the cluster centers). See [l^l for a discussion of 
redshift space effects interpreted using a homogeneous 
model, especially Fig. 5 illustrating caustics related to the 
turnaround radius. Velocity dispersions of 91 rich clus- 
ters of the ESO Nearby Abell Cluster Survey (ENACS) 
and earlier surveys selected out to z = 0.1 from a solid 
angle of w 8400 deg^ have a mean of (t„ ~ 642 ± 24 km/s 
(Table la, [i^; standard error in the mean; cf Fig. 5a of 
jis'l). An SDSS Data Release 4 analysis of 109 clusters 
out to z < 0.1 over 6700 deg^ yields = 585 ± 15 km/s 
(Table 1, l49|; standard error in the mean). The two cat- 
alogues have very few members in common, so an inverse- 
variance weighted mean can be calculated: 



ay = 601 ± 13km/s. 



(36) 



Thus, with ft. = 0.74 from 
(l36t we set 



TTCom 

pec 



(0) := 



2n 



infall 



and using ([M)) . (I55t . and 

4cr.„ 



^void ^void 

= 36±3km/s/Mpc. 



(37) 



D. Effective metric 



•S E 

a 



FIG. 2. Effective zero redshift matter density parameter 
finf (|32p in the virialization approximation. The upper, cen- 
tral, and lower curves correspond to -ffpec"(0) = 33,36, and 
39 km/s/Mpc, respectively [see p7l) ]. 



and c is the conversion constant between space and time 
units. Equation (j39p can be integrated numerically to 
obtain and the luminosity distance 



z)Rfsinh 



Dcff ■ 



(40) 



The effective differential volume element per square de- 
gree, necessary for faint galaxy number count analyses, 
is 



dF 



eff 



dz 



f— V 

Vl80/ 



eff 



dx' 



eff 



peff 



dz 



(41) 



The luminosity distance for the virialization approxi- 
mation and for the homogeneous EdS and ACDM models 
can be normalized to the Milne model for convenience, 
yielding distance moduli m — M. The fraction of the dis- 
tance modulus that would be needed to correct the ho- 
mogeneous EdS model to match the ACDM model can 
be written as follows: 



/r 



logipdf -logipdfs 
logio4CDM_iog^^rfEdS- 



(42) 



We extend the usual FLRW metric to an effective, 
spherically symmetric, comoving, observer-centred ex- 
pression for the metric: 



ds' 



eff^ 



sinh' 



eff 



where the differential radial comoving distance is 
cdt 



cdt 



a da 



■ da 



a- 



H°«{z) 



da. 



(38) 



(39) 



III. RESULTS 

The virialization fraction /vir from the two iV-body 
simulations has already been shown in Fig. [TJ Figure [2] 
shows that, as expected, the effective matter density pa- 
rameter contributing to curvature drops as the redshift 
decreases to zero. The low redshift value, rj^f (0) = 
0.26 ±0.05, is impressively close to observational low red- 
shift estimates of ilm ~ 0.3. It is not a result of fitting the 
virialization approximation to observational data. Apart 
from assuming an EdS background model, the only ob- 
servational values used in the formulae above are those 



presented in Sect. [nCl i.e. H°^{{)) ^ and 



rrcom 
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z 

fvir 



TABLE I. A''-body simulation estimates of fvir{z) (U) 



10.0 5.0 3.0 2.0 1.5 1.0 0.5 0.3 0.1 0.0 

1.4 X 10"^ 0.0036 0.026 0.061 0.12 0.22 0.27 0.32 0.35 




FIG. 3. 



Ratio of effective to background expansion rates 
using (I24|, ^Jii, and (HH). The upper, cen- 
tral, and lower curves correspond to ffp°™(0) = 39,36, and 
33 km/s/Mpc, respectively. 




FIG. 4. Effective comoving curvature radius Rq (I3ip in the 
virialization approximation. The uncertainty in /fp°™(0) is 
not visible in this plot. 



([57)1 . Figure [3] shows the effective expansion rate evolu- 
tion, which, along with /vir, explains the lower effective 
matter density parameter at low redshifts. Through the 
comoving curvature radius ()3ip . the void-dominated ex- 
pansion rate at low redshifts also leads to stronger neg- 
ative average curvature, by decreasing Kq . Figure [4] 
shows the effective curvature radius evolution R'q[z). 

Figure [5] shows that applying the virialization approxi- 
mation to the EdS model brings the EdS distance modu- 
lus close to the (homogeneous) ACDM distance modulus. 
The fraction fm~M (|42p provided by the virialization ap- 
proximation is shown in Fig. [6l showing that a discrep- 
ancy remains for z < I. Figure [7] shows that the volume 



o 




FIG. 5. (colour online) Distance moduli for the EdS 
(bottom curve) and ACDM (top curve for z < 0.5, black 
online) homogeneous models, and for the EdS virializa- 
tion approximation (Sect. [TTl green online, thick curve for 
HpT{0) = 36 km/s/Mpc, thin curves for H;°T{0) = 33 and 
39 km/s/Mpc), all normalized to the Milne model (JlmO = 
O^^Ao). 




FIG. 6. Fraction of EdS-to-ACDM distance modulus fm~M 
(|42p provided by the virialization approximation. 



element needed to match faint number counts is provided 
by this approximation. 
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FIG. 7. (colour online) Differential volume element dV" /dz 
(gg for the EdS (bottom curve) and ACDM (black online) 
homogeneous models and for the EdS virialization approx- 
imation (thick and very thin curves for Hp°^{0) = 39,36, 
and 33 km/s/Mpc, from top to bottom, respectively, green 
online), in comoving Mpc"^/sq.deg. 



IV. DISCUSSION 

The dominance of negative curvature over positive cur- 
vature on large scales at recent epochs has already been 
postulated based on detailed kinematical and curvature 
backreaction estimates [11] and models 
fssj . By definition, over-and under-densities of comoving 
perturbations in an FLRW background model average 
to zero. In general, this is also true on mass-preserving 
Lagrangian domains V in the nonlinear regime. The 
restmass conservation law assures compensation to hold 
throughout the evolution of structure. However, the in- 
trinsic curvature does not obey a conservation law: only 
a dynamical relation between the kinematical backreac- 
tion, curvature, and volume evolution is conserved p3| . 
as is used in this paper through ^ with Ha = 0; see the 
discussion in (s^ . 

Virialized regions are, by definition, highly concen- 
trated, occupying very little volume. Thus, while the 
average of the density deviations from the background 
tends to zero, the average peculiar curvature (average 
curvature with respect to the background FLRW model 
curvature) tends to a negative value when a high fraction 
of matter has virialized. The key formula in [23| is (104), 
with a general discussion in Section 5 of strategies for 
calculating a more detailed approximation than the one 
presented here. 

At small enough scales, the backreaction model based 
on a volume average of the relativistic Zel'dovich approx- 
imation mentioned above, representing the effects of the 
statistical formation of structure in detail [24, ,25], shows 
that positive curvature dominates in collapsing regions. 
It is speculated that, on these smaller scales, positive cur- 
vature plays the role that is usually attributed to "dark 
matter" , i.e. that a substantial portion of dark matter 
might also be an artefact of using a Newtonian approxi- 



mation to structure formation. Although the estimates in 
this paper deal with void scales and larger, detailed rela- 
tivistic calculations would have to take this into account. 
We do not attempt to model this here, since we consider 
our results to be an initial quantitative estimate of an im- 
portant physical effect, suggesting that it is worthwhile 
to go beyond the model presented here. 

For over a decade, radially inhomogeneous models, 
both Stephani 56] and LTB [53, HI] models, have been 
proposed as dark-energy-free, relativistic fits to the type 
la supernovae luminosity-distance-redshift relation, and 
have become best known as "void models" [59|, [GO] on 
a scale of up to ~ 2.5 Gpc (e.g., although "hump 

models" have also been inferred from the data [2^, [s^ . 
There has been much debate over whether or not the fits 
should be excluded on either observational grounds or as 
bein g in conflict with the Cosmological Principle. (See, 
e.g., 11.) 

The virialization approximation (which is neither a 
Stephani nor an LTB model) suggests another interpre- 
tation that may resolve some issues in these debates. It is 
true that in a comoving, constant time slice at the present 
epoch to, SL reasonable projection of our past light cone 
observations forward to to is hard to reconcile — at least 
for large-scale homogeneous models — with our location 
near the centre of a void as large as a gigaparsec. How- 
ever, cosmological observations are not made on the to 
time slice. On the past light cone, the virialization epoch 
approximately corresponds to an observer-centred "void" 
of a gigaparsec or so in size, where "void" here means, 
e.g., the spherical region around the observer extending 
over the redshift range where fvh-{z) ^0.1. The "void" 
that is relevant for an effective metric is that defined by 
the volume-averaged density parameter on our past light 
cone, and this "void" is mirrored by virialization. We 
are located at a highly privileged, highly non-random, 
centralized position on the past light cone: at the apex 
of the cone, which is the epoch of highest virialization. 
Thinking spatially at constant t, locating the observer 
at the centre of a (non-averaged) void is a hypothesis, 
while on the past light cone, our privileged position is a 
geometrical fact. 

With hindsight, the numerical results of LTB void 
models by, e.g., (63{ and the LTB void model Ho{r) es- 
timates of [6J] , for a radial coordinate r at constant cos- 
mological time, roughly correspond to the results found 
here, but with fundamental differences in their deriva- 
tion and interpretation. Here, we start with a large- 
scale homogeneous (EdS) model and use standard A^- 
body simulations to quantify an approximate, effective 
metric at low redshifts, providing a general-relativistic 
correction to the FLRW metric. Our resulting "void" 
is a past-light-cone, virialization-epoch, volume-averaged 
pseudo-void only. In the LTB models, Hq and J7niO vary 
with the radial coordinate r at a constant cosmologi- 
cal time t, at which there is a true void rather than a 
virialization-epoch "void" . However, the LTB model pa- 
rameters are motivated by and interpreted in relation to 
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FIG. 8. Effective deceleration parameter q"^ (1431) for the 
EdS virialization approximation presented in this paper. 

observations — on the past light cone. Thus, it is unsur- 
prising that the numerical results are roughly similar to 
ours. 

A more closely related inhomogeneous, dark-energy- 
free relativistic model is the timescape model (formerly 
known as the "fractal bubble" model) [65l - [67j . At 
late epochs, this is a negatively-curved, void-dominated, 
multi-scale, scalar averaging model that focuses particu- 
larly on time calibration. The model compares well to the 
ACDM model in fitting supernovae type la data [6^ and 
has been used to interpret bulk velocity flows on scales 
up to 65 Mpc from the Local Group [69j . 

Is the scale factor accelerating according to the viri- 
alization approximation? The study of LTB models 
shows that for inhomogeneous universe models in general, 
at least two different definitions of the deceleration pa- 
rameter can be usefully made, and luminosity-distance- 
redshift observations do not imply model-free estimates 
of either of these [70!|. Here, using the path of a radial 
photon in (pS)) . dx°^/da from and appropriately 

converting between space and time units, an effective de- 
celeration parameter. 



can be defined. Figure |S] shows that for this definition, 
moderate acceleration occurs for z ^ 1.5. 

Given that only two observational values are used in 
the above approximation, i.e. one that is commonly de- 
rived from observations, iJ°^(0) (1^51) . and one that we 
derive here from observations, iJp°™(0) (|371) . it would 
be to see if the approximation agreed with other further 
observational constraints beyond f^J^(O) = 0.26 ± 0.05 
and the d'^^{z) relation. One important constraint is the 
age of the Universe at the (approximately comoving) ob- 
server's spacetime location. Using the EdS age at z = 10, 
the virialization approximation gives the present-day age 
of the Universe to = 12.2tp;4 Gyr. This is a little low, 
but consistent with to = 13^2 Gyr from the stellar popu- 
lation dating of moderate redshift elliptical galaxies over 



0.3 < z < 0.9 OT to> 12.8 ± 1 Gyr from globular 
cluster ages calibrated using Hipparcos parallaxes (72i] . 

V. CONCLUSION 

The approximation presented above assumes "old 
physics" (general relativity), not "new physics". It is 
a rough approximation of old physics, general relativ- 
ity, applied to the old observations of a high virialization 
fraction (galaxies, clusters of galaxies) at recent epochs 
and at a small distance scale in the past light cone. Nu- 
merical values of /vir(-z) are esimated from a standard 
EdS A^-body simulation [l^, [l^ . Only two observational 
values are assumed: the zero-redshift Hubble parameter 
(PS]) and the zero-redshift peculiar expansion rate across 
non-virialized regions, i?p°™(0) ([57]) . The inferred effec- 
tive low redshift matter density parameter is realistic, 
^nf = 0.26 ± 0.05, and the virialization-corrected EdS 
distance modulus is close to the ACDM distance modu- 
lus (Figs m El). 

Of course, virialized objects cannot be completely ig- 
nored, the virialization fraction is derived from an A^- 
body simulations with a fixed (zero) curvature EdS 
model, and different choices of parameters in the A^- 
body analysis (group finder algorithm, group finder min- 
imum number of particles, definition of the void size) and 
the A^-body simulation itself (calculation algorithm, vol- 
ume of the simulation, particle mass resolution) would be 
likely to modify the above calculations. Our proposed in- 
terpolation (|^ and normalization ([57)1 of ifp°™(z) could 
also be improved in many ways. It is, however, un- 
likely that these refinements would substantially reduce 
the amplitude of the corrections estimated here, since 
the observational evidence in favor of a high virialization 
fraction at the present epoch z <^ 1 and the observational 
and theoretical evidence for H{z) + H^°'^{z) a'^ ^ to 
roughly follow the values shown in Fig. [3] is strong. A 
more detailed approximation than that presented here 
might show that a hyperbolic or spherical, zero-dark- 
energy background FLRW metric, corrected for virializa- 
tion, would better match the full range of observational 
constraints. It may also be important to consider inho- 
mogeneous light-path effects on standard candles, such as 
type la supernovae, as a bias when assuming the FLRW 
models as a family of background models |73l - [75| . High 
redshift observations should be modeled and interpreted 
in a way that avoids dependence on the low redshift 
metric. Low-redshift peculiar-velocity observations (e.g. 
fld^ ) would be useful if designed to optimally estimate 

Pending more accurate, relativistic calculations, it 
would seem prudent to consider "dark energy" as an 
artefact of virialization-induced negative spatial curva- 
ture and void-dominated expansion rates, both of these 
being physical properties that are neglected in the stan- 
dard cosmological framework. 



10 



ACKNOWLEDGMENTS 

Thank you to Martin Kerscher, Tomasz Kazimierczak 
Alexander Wiegand, Bruce A. Peterson, and Bartosz 
Lew for useful comments. Some of this work was car- 
ried out within the framework of the European Asso- 
ciated Laboratory "Astrophysics Poland-France". BFR 
thanks the Centre de Recherche Astrophysique de Lyon 
for a warm welcome and scientifically productive hos- 
pitality. A part of this work was conducted within 
the "Lyon Institute of Origins" under grant ANR-10- 
LABX-66. Some of JJO's contributions to this work 
were supported by the Polish Ministry of Science and 
Higher Education under "Mobilnosc Plus II edycja" . A 



part of this project has made use of Program Obliczeh 
Wielkich Wyzwan nauki i techniki (POWIEW) compu- 
tational resources (grant 87) at the Poznah Supercom- 
puting and Networking Center (PCSS). The TV-body 
simulations analyzed in this paper were carried out by 
the Virgo Supercomputing Consortium using computers 
based at the Computing Centre of the Max-Planck So- 
ciety in Garching and at the Edinburgh Parallel Com- 
puting Centre [77|. Use was made of the GNU Oc- 
tave command-line, high-level numerical computation 
software (http://www.gnu.org/software/octave), and 
the Centre de Donnees astronomiques de Strasbourg 
( jhttp : / / cdsads . u- strasbg . f rj . 



M. Fukugita, K. Yamashita, F. Takahara, and Y. Yoshii, 

ApJL 361, LI (Sep. 1990) 

Y. Yoshii and B. A. Peterson, |Astrophys. J.| 444, 15 (May 
1995) 

M. Chiba and Y. Yoshii, Astrophys. J. 489, 485 (Nov. 
1997) 

B. Fort, Y. Mellier, and M. Dantel-Fort, A&A 321, 353 
(May 1997) 

S. Perlmutter, G. Aldering, G. Goldhaber, et al., Astro- 
phys. J. 517, 565 (Jun. 1999), arXiy:astro-ph/98 12133 
B. Suntzeff, M. M. Phillips, 
46 (Nov. 1998), 



P. Schmidt^ 
al.. 



N. 



507, 



ApJSupp 



Astrophys. J. 
[arXiv:astro-ph/9805200 

D. N. Spergel, L. Verde, H. V. Peiris, et al., 
148, 175 (Sep. 2003), astro-ph/0302209 

Planck Collaboration, P. A. R! Ade, N. Aghanim, 
C. Armitage-Caplan, M. Arnaud, M. Ashdown, F. Atrio- 
Barandela, J. Aumont, C. Baccigalupi, A. J. Banday, and 
et al., ArXiv e-prints(Mar. 2013), arXi v:1303. 50761 
W. de Sitter, MNRAS 78, 3 (Nov. 1917) 
A. Friedmann, Mir kak prostranstvo i vremya (The Uni- 
verse as Space and Ttme) (Petrograd: Academia, 1923) 
A. Friedmann, Zeitschr. fiir Phys. 21, 326 (1924) 

G. Lemaitre, Annales de la Societe Scientifique de Brux- 
elles 47, 49 (1927) 

H. P. Robertson, Astrophys^ J. 82, 284 (Nov. 1935) 

T. Buchert, General Relativity and Gravitation 40, 467 
(Feb. 2008), ar Xiv:0707.2153l 

E. W. Kolb, V. Marra, a nd S. Matarrese, 
[General Relativity and Gravitation| 42, 1399 (Jun. 
2010), "arXiv:0901.4566 

J. Larena, J.-M. Alimi, T. Buchert, M. Kunz, and P.- 
S. Corasaniti, Phy's. Rev7D, 79, 083011 (Apr. 2009), 
arXiv:0 808. 116 11 

T . Buchert, G. F. R. Ellis, and H. van Elst, 
.General Relativity and Gravitationl 41, 2017 (Sep. 2009), 
[arXiv:0906.0134 

A. Jenkins, C. S. Frenk, F. R. Pearce, P. A. Thomas, 
J. M. Colberg, S. D. M. White, H. M. P. Couch man, J. A . 
Peacock, G. Efstathiou, and A. H. Nelson, Astrophys. J.| 
499, 20 (May 1998), arXiv:astro-ph/9709010T 
P. A. Thomas, J. M. Colberg, H. M. P. Couchman, et al., 
[MNRAS 296, 1061 (Jun. 1998) 

J. P. Ostriker and P. J. Steinhardt, ArXiv Astrophysics 



e-prints(May 1995), |arXiv:astro-ph/9505066| 
[21] C. Lacey and S. Cole, MNRAS 262, 627 (Jun. 1993) 
[22] T. Buchert, M. Kers cher, and C. Sicka, Phys. Rev. D 62, 

043525 (Aug. 2000), [arXl7:astro-ph/9912347| 
[23] T. Buchert and M. Carfora, 

I Classical an d Quantum Gravity] 25, 195001 (Oct. 

2l508)7axXiv:0803.1401 

[24] T. Buchert and M. Ostermann, Phys. Rev. D| 86, 023520 

(Jul. 2012), arXiv: 1203.6263 
[25] T. Buchert, C. Na yet, and A. W iegand, Phys. Rev. D, in 

press (Mar. 2013). [ajXrv:1303.6 193 
[26] G. Lemaitre, Annales de la Societe Scientifique de Brux- 

elles 53, 51 (1933) 

[27] R. C. Tolman, Proceedings of the National Academy of Science| 

20, 169 (Mar. 1934) 
[28] H. Bondi, MNRAS 107, 410 (1947) 

[29] M. Celerier, K. Bolejko, and A. Krasinski, 'A&A 518, 

A21 (Jul. 2010), arXiv: 0906. 0905 
[30] E. W. Kolb and C. R. Lamb, ArXiv e-prints(Nov. 2009), 

arXiv:0911.3852 
[31] |http : / /arXiv . org/abs/1303 . 4444 

[32] B. F. Roukema, P. J. Quinn, and B. A. Peterson, 
in Observational Cosmology: an International Sympo- 
sium, Milano, Italy, 21-25 September 1992, Astronom- 
ical Society of the Pacific Conference Series, Vol. 51, 
edited by G. L. Chincarini, A. lovino, T. Maccacaro, & 
D. Maccagni (1993) p. 298 

[33] B. F. Roukema, B. A. Peterson, P. J. Quinn, and 
B. Rocca-V olmerange, MNRAS 292, 835 (Dec. 1997), 
|a rXiv:astro -ph/9707294 

[34] C. Lacey and S. Cole, in Observational Cosmology: an 
International Symposium, Milano, Italy, 21-25 Septem- 
ber 1992, Astronomical Society of the Pacific Conference 
Series, Vol. 51, edited by G. L. Chincarini, A. lovino, 
T. Maccacaro, and D. Maccagni (1993) p. 192 

[35] B. F. Roukema, S. Ninin, J. Devriendt, F. R. Bouchet, 
B. Guiderdoni, and G. A. Mamon, A&A 373, 494 (Jul. 
2001), arXiv:astro-ph/0105152 

[36] http : / /www-hpcc . astro . Was hington. edu/toolsT] 
f of .html] 

[37] J. Tinker, A. V. Kravtsov, A. Klypin, K. Abazajian, 
M. Warren, G. Yepes, S. Gottlob er, and D . E. Holz, 
Astrophys. J. 68 8, 709 (Dec. 2008), arXiv:0803. 2706 

[38] M. White, lAfcAl 367, 27 (Feb. 2001), 



11 



143, 241 (Dec. 2002), 



arXiv:astro-ph/0011495" 
M. White, ApJSupp 
|arXiv:astro-ph/ 0207185 
A. Wiegan d and T. Buchert, |Phys. Rev. D 82, 023523 
(Jul. 2010). [arXivT 002.3912 

P. J. E. Peebles, Large-Scale Structure of the Universe 
(Princeton University Press, — cl980, 1980) 
A. G. Riess, L. Maori, S. Casertano, H. Lampeitl, 
H. C. Ferguson, A. V. Filippenko, S. W. Jha, W . Li, 
R. Chornock, and J. M. Silverman, As trophys. J.| 730, 
119 (Apr. 2011), arXnFril03. 2976 

W. L. Freedman, B. F. Madore, V. Scowcroft, C. Burns, 
A. Monson, S. E. Persson, M. Seibert, and J. Rigby, 
Astrophys. J. 758, 24 (Oct. 20 12), arXiv :1208. 3281 
M. Pl ionis and S. Basilakos, JMNRASI 330, 399 (Feb. 
2002), 'arXiv:astro-ph/0106491 



C. Foster and L. A. Nelson, Ast rophys. TT] 699, 1252 (Jul. 
2009), arXiv:0904.4721 ^ 

K. Rines, M. J. Geller, M. J. Kurtz, and A . Diaf erio, .AJ. 

126, 2152 (Nov. 2003), arXiv:astro-ph /0306538] 

N. Kaiser, MNRAS 227, 1 (Jul. 1987) 

A. Mazure, P. Katgert, R. den Ha rtog, et al., A&A 310, 

31 (Jun. 1996), arXiv:astro-ph/9511052 

J. A. L. Aguerri, R. Sanchez- J anssen, and C. Mufioz- 

Tuiion, A&A 471, 17 (Aug. 2007), arXiv:0704.1579 

S. Rasanen, Classical and Quantum Gravity 28, 164008 

(Aug. 2011), arXiv: 1102.0408 

T. Buchert, Classical and Quantum Gravity 28, 164007 
(Aug. 2011), arXiv:1103.2016 

T. Buchert and S. Rasanen, 
[Annual Review of Nuclear a nd Particle Sciencel 62, 
57 (Nov. 2012), arXiv:1112.53^ 

D. L. Wiltshire, Classical and Quan tum Gravity| 28, 

164006 (Aug. 2011), arXiv:1106.1693 

X. Roy and T. Buchert, Classical and Quan tum Gravity] 
27, 175013 (Sep. 2010), arXiv:0909.4155 

C. Clarkson, G. Ellis, J. Larena, and O. Umeh, ArXiv 

e-prints(Sep. 2011), arXiv:1109.2314 

M. P. Dabrowski a nd M. A. H endry, A strophys. J.| 498, 

67 (May 1998), arXiv:as tro-ph/9704123| 

N. Mustapha, C. HeUab y, and G. F. R. E Uis, MNRAS 

292, 817 (Dec. 1997), arXiv:gr-qc/9808079 

M.-N. Ce lerier, ArXiv A strophysics e-prints(Jul. 1999), 

|arXiv:astro-ph/9907206] 



[59 

[6o; 

[61 

[62; 

[63 

[64; 
[65; 

[66 
[67 
[68; 
[69 

[7o; 

[71 

[72; 
[73; 

[74; 

[75; 

[76; 

[77 



K. Tomita, MNRAS 326, 287 (Sep. 2001), 
arXiv:astro-ph/00114 84 

K. Bolejko, M.-N. Celerier, and A. Krasihski, 
'Classical and Quantum G ravity, 28, 164002 (Aug. 
2011), arXiv:1102.1449 

J. Garcia-BeUido and T. Haugb0lle, [jClPj 4, 003 (Apr. 
2008), arXiv:0802.1523 

A. Krasinski and K. Bolejko, ArXiv e-prints(Dec. 2012), 

arXiv: 12 12,4697 

H. Alnes, M. Amarzgui oui, and 0. Gr0n, |Phys 



Rev. D| 



451 



73, 083519 (Apr. 2006), arXiv:astro-ph/0512006 
K. Enqvist, General Relativity and Gravitation 40 

(Feb. 2008), arXiv:0709.20 44j 

D. L. Wihshire, New Jo urnal of Physics] 9, 377 (Oct 



2007), arXiv:gr-qc/0702082j 

D. L. Wihshire, Physical Re view Letters] 99, 251101 
(Dec. 2007), arXiv:0709.0732 

D. L. W ihshire, Phys. Rev. D 80, 123512 (Dec. 2009), 
larXiv:0909.0749 

P. R. Smale and D. L. Wihshire. IMNRAS1 413. 367 (May 
2011), arXiv:1009.5855 

D. L. Wiltshire, P. R. Smale, T. Mattsson, and 
R. Watkins, ArXiv e-prints(Jan. 2012), arXiv:1201.5371J 
K. Bolejko and L. Andersson, JCAP 10, 003 (Oct. 2008), 
arXiv:0807.3577 

I. Ferreras, A. Me lchiorri, and J. Silk, MNRAS 327, L47 

(Nov. 2001), arX iv:astro-ph/01 05384 

L. M. Kr auss, Ph ys RepTr 333, 33 (Aug. 2000), 

]arXiv:ast ro-ph/9907308] 

K. Kainulainen and V. Marra, in 
\Ame rican Institute of Physics Conference Series , Amer- 
ican Institute of Physics Conference Series, Vol. 1241, 
edited by J.-M. Alimi an d A. Fuozfa (2010) pp. 
1043-1050, larXiv:0911.5"584l 

C. Clarkson, G. F. R. Ellis, A. Faltenbacher, 
R. Maartens, O. Umeh, and J.-P. Uzan, JMNRAS 426, 
1121 (Oct. 2012), arXiv:110a2484 

P. Fleury, H. Dupuy, and J.-P. Uzan, ArXiv e-prints(Feb. 
2013) , larXiVnM23308l 

R. B. TuUy, E. J. Shaya, L D. Karachentse v, H. M. Cour- 
tois, D. D. Kocevski, L. Rizzi, and A. Peel, ]Astrophys."jT] 
676, 184 (Mar. 2008), arXiv:0705.4139 
http : //www ■ mpa-garching . mpg . de/NiunCos ] 



